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Experiments on Bose-Einstein condensation (BEC) of excitons in the semiconductor CU2O started 
over thirty years ago, as one of the first serious attempts at exciton BEC. Early claims were based 
on spectroscopic signatures and transport data which have since been reinterpreted, in large part 
because the Auger recombination process for excitons was not well understood. Understanding 
of the Auger process has also advanced, and recent experiments have made significant progress 
toward exciton BEC. We review the history of experiments on exciton BEC in CU2O, the Auger 
recombination process, and the prospects for observing exciton BEC in this system in the near 
future. 

PACS numbers: 



I. INTRODUCTION 

Following the initial predictions of Bose-Einstein con- 
densation (BEC) of excitons by Moskalenko 1] , Blatt and 
coworkers 0, and Casella Q in the 1960's, and founda- 
tional theory by Keldysh and Kozlov [3, [EI and Hanamura 
and Haug [6J-18JI. (with significant later work by Comte 
and Nozieres [3, [l(| and others [ll|), experimental re- 
search on Bose-Einstein condensation of excitons began 
in the 1970's with a focus on deeply bound excitons in 
bulk semiconductors. Deeply bound excitons were espe- 
cially attractive because they held out the possibility of 
a room temperature condensate. The critical tempera- 
ture for BEC increases with increasing density, so high 
exciton densities are favored, but when the exciton den- 
sity n approaches the critical value n ~ a _3 , where a is 
the exciton radius, the Fermi nature of the constituent 
electrons and holes in the excitons becomes important, 
and a transition to a degenerate electron-hole plasma can 
occur (see Ref. Section 2.5.4.) Since the radius of 
the excitons is smaller when the excitons are more tightly 
bound, materials with low dielectric constant and tightly 
bound excitons are favored. A large exciton binding en- 
ergy is also favored for high temperature studies because 
excitons can disassociate into free electrons and holes at 
temperatures comparable to their Rydberg energy (l3| . 

The earliest experiments on excitonic condensation in 
bulk semiconductors were done with the bulk semicon- 
ductors CdSe and CuCl The most complete 
study on CuCl was done in the early 1980's [H, [ljj. In 
CuCl, there is a tightly bound biexciton state (i.e., an 
excitonic molecule like H2), which, of course, is also an 
integer-spin boson which can in principle undergo BEC. 
The biexciton binding energy in CuCl is 26 meV, which 
is greater than fc^T at room temperature. The experi- 
ments with CuCl ran into interpretive difficulties, how- 
ever, because CuCl has a strong polariton effect, in which 
photons in the medium couple strongly to the exciton or 



biexciton states, leading to a very short radiative lifetime 
for low-momentum excitons. At high density, the biexci- 
tons in CuCl showed a distinct peak in low-energy states, 
but were far from equilibrium. The enhancement of the 
population in low-energy states was probably related to 
the boson nature of the particles through the process of 
stimulated scatte ring, w hich drives the onset of BEC in 
all bosonic gase s I2li 21 1. In recent years, several sophis- 
ticated models [22H24j | have been developed which de- 
scribe the nonequilibrium behavior of polaritonic gases 
and show clearly the role of the bosonic stimulated scat- 
tering. It would be quite interesting to apply these mod- 
ern theoretical tools to the early CuCl experiments to 
estimate to what degree the occupation number of the 
biexcitons was enhanced by their bosonic nature. 

Around the same time, the semiconductor CU2O was 
identified as a good candidate for excitonic BEC. The ex- 
citons in CU2O have binding energy of 150 meV [25|]; the 
biexciton state is either nonexistent or only very weakly 
bound, on the order of 1 meV or less [2y, [27]], so that 
one can generally ignore the existence of biexcitons (we 
discuss this further below). Furthermore, the lowest ex- 
citon state in CU2O, known as the paraexciton, has no 
polariton effect, because it has zero oscillator strength 
for interaction with photons, due to the selection rules 
for this state which arise from the symmetry of the crys- 
tal [25| . In principle, this means that these excitons have 
infinite lifetime, but a weak phonon-assisted recombina- 
tion process is allowed [28], and also recombination at 
impurities can occur, so that the lifetime of paraexcitons 
in Cu 2 is typically of the order of hundreds of nanosec- 
onds to microseconds [29l - [3l| . This is very long com- 
pared to excitons in most semiconductors, and very long 
compared to the thermalization times of the excitons via 
exciton-phonon interaction, which is tens of picoseconds 
[32j. Because there has never been a large technological 
effort to produce pure CU2O, most of the purest samples 
come from rock collectors who obtain them from copper 
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mines. 

It is important to note the difference of exciton BEC 
(EBEC), from a very different state known as the 
electron- hole liquid (EHL), sometimes also called "con- 
densation." This state received a significant amount of 
attention in the 1970's and 1980's 0, EH following its 
observation in two bulk semiconductors, Ge and Si. This 
state is analogous to liquid mercury: the electrons and 
holes generated by optical excitation are not bound into 
pairs as excitons or biexcitons but instead form two inter- 
penetrating Fermi gases with the properties of a classical 
liquid, with a surface tension. The EHL state is a con- 
ductor, while the EBEC state is an insulator, because 
excitons are charge neutral. 

The main difference between semiconductors like Ge 
and Si, which allow EHL formation, and a semiconductor 
like Cu 2 0, which does not allow EHL but allows EBEC, 
is the fact that Si and Ge are indirect-gap semiconduc- 
tors with multiple degenerate valleys in the conduction 
band. This degeneracy allows the density of the excited 
electrons to be high while still keeping the average ki- 
netic energy of the carriers low. This means that the 
mutual attraction of the electrons and holes overcomes 
the kinetic energy cost of the Fermi level of free carriers, 
leading to a net free energy savings to enter the EHL state 
at high excited carrier density |34j . In a nice experiment, 
Timofeev and coworkers [35[ used a stress geometry to 
lift this degeneracy in a germanium crystal, and magnetic 
field to prevent biexciton formation, and found that the 
EHL state did not occur, leaving only free excitons. At 
high density they saw evidence for Bose- Einstein statis- 
tics of the excitons (see Ref. [ll|, section 1.4.2), although 
not for EBEC. 

A very important step was accomplished in the 1980's 
in the trapping of excitons in CU2O. Inhomogeneous 
stress was used to create a harmonic potential minimum 
for the excitons inside the crystal [3a . very similar 
in geometry to the magneto-optical traps used to hold 
cold atoms for BEC experiments. As discussed below, 
some early experiments on EBEC used surface excita- 
tion of the crystal, which had the advantage of achieving 
high carrier density due to the short absorption length of 
the excitation light, but generated a highly nonequilib- 
rium situation with an expanding gas that was difficult 
to analyze. The stress used in the exciton traps has the 
side effect that it changes the symmetry of the crystal, 
so that the paraexciton state no longer has a forbidden 
radiative recombination transition [38j, but the lifetime 
of the paraexcitons in a stressed crystal is still very long 
compared to the exciton-phonon interaction time. 

The phonon-assisted recombination process in CU2O 
provides nice data to indicate the behavior of the exci- 
tons and whether they undergo BEC. This process can 
occur for both the paraexciton, which is a singlet state, 
and the orthoexciton state, which is a triplet that lies 12 
meV above the paraexciton in unstressed crystals. The 
orthoexciton state also has a symmetry-allowed radiative 
recombination process. (The singlet-triplet terminology 
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FIG. 1: Solid lines: luminescence from excitons in CU2O from 
low temperature to room temperature. Dashed line: fit to 
the theory for phonon-assisted luminescence discussed in the 
text, giving the thermal distribution of the excitons. A single 
parameter, namely the lattice temperature, is used to fit all 
the curves. From Ref. [53] ■ 



of CU2O sometimes confuses those who work with or- 
ganic light emitters, in which the state with forbidden 
recombination is called the triplet state and the state 
with allowed radiative recombination is called the singlet 
state. The difference in the terminology is because in 
the case of organics, the transition is assumed to have an 
s-p dipole-allowed matrix element, and the spin singlet 
state corresponds to no spin flip of the excited electron, 
and thus an allowed process, while the spin triplet state 
corresponds to a spin flip. In the case of CU2O, both 
the conduction and valence bands are s-like, so that the 
transition is forbidden for the paraexciton, but the spin 
triplet state makes the transition allowed by symmetry.) 

Figure [T] shows typical spectra of phonon-assisted lu- 
minescence from orthoexcitons in CU2O at various tem- 
peratures up to room temperature. Because an optical 
phonon is emitted or absorbed when an exciton recom- 
bines, an exciton at any momentum can recombine via 
this process, with the phonon taking up any excess mo- 
mentum. This is in contrast to the standard direct re- 
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combination process for excitons in most semiconductors, 
such as GaAs, which is allowed only for excitons with 
low momentum — the exciton momentum and the emitted 
photon momentum must match in the direction parallel 
to the surface of the crystal. For the phonon-assisted 
process in CU2O, since the energy of the optical phonons 
is nearly constant, the energy of the emitted photon is 
equal to the total energy of the exciton before it recom- 
bined minus a constant energy, namely the energy of the 
optical phonon. The energy spectrum of the phonon- 
assisted luminescence therefore gives the kinetic energy 
distribution of the excitons directly. If we take the matrix 
element as independent of the exciton momentum, which 
is correct in CU2O to a high degree, the intensity of the 
light emitted at a given energy is directly proportional 
to the number of excitons with the corresponding kinetic 
energy. The theory fits to the data in Fig. [1] are therefore 
given simply by I(E + E ) oc D(E)f(E), where E is the 
exciton kinetic energy, D(E) is the density of states of the 
excitons, proportional to E l l 2 in three dimensions, and 
f(E) is the occupation number of the excitons, propor- 
tional to e - E / k sT a ^ i ow density There are two phonon- 
assisted lines with the same shape, one for emission of 
an optical phonon, and one for absorption of an optical 
phonon. The relative height of these two lines depends 
on the temperature of the lattice. At low-temperature, 
the phonon-emission line dominates, since there are no 
phonons to absorb, while at high temperature, the two 
lines become comparable. 

When many-body interactions are taken into ac- 
count, the analysis of the phonon-assisted luminescence is 
slightly more complicated. Each momentum state emits 
with an energy profile given by the spectral function 
A(k,u)), which is approximately equal to 5(ljJ — uit) at 
low density and low temperature, but can be substan- 
tially broadened when there are strong particle- particle 
interactions [H, IM, E3| • This brings up a crucial differ- 
ence between exciton BEC experiments and experiments 
on cold atoms. The fact that excitons can recombine and 
turn into photons allows us to have direct access to the 
spectral function. (This property has been used widely 
in the exciton-polariton BEC experiments,_e.g. to show 
spectral narrowing due to coherence [4l|, S3-) In c °ld 
atom experiments, the spectral function is a theoretical 
construct that is never directly observed. 

Direct access to the spectral function means that we 
can see absolute shifts of the ground-state energy of the 
excitons directly as energy shifts of the emitted photons, 
and that we can see the effect of broadening of the spec- 
tral function as spectral broadening of the photon emis- 
sion. To lowest order, the former can be identified as the 
real part of the exciton self-energy, and the latter as the 
imaginary part of the self-energy. (See Ref. [T2], section 
8.4). In Fig. [TJ the exciton-phonon interaction increases 
as the temperature increases, leading to a red shift of the 
luminescence line due to renormalization of the real self- 
energy and broadening due to increasing imaginary self 
energy. The imaginary self-energy is, to first order, just 



given by T = h/r, where r is the average scattering time 
of the excitons. 

This effect was used in the first attempt to estimate 
the absolute density of the excitons independently from 
the assumed statistics of the excitons in the CU2O ex- 
periments [43j . as discussed below. Given an estimate of 
the exciton-exciton scattering cross section, the scatter- 
ing time r for a given density can be calculated for the 
case when the lattice is cold, near 1 K, so that exciton- 
phonon interactions are negligible, and this can be com- 
pared to the observed broadening as a function of exciton 
density. 

The cross section for exciton-exciton scattering is ac- 
tually a problematic theoretical topic for exciton physics, 
as compared to cold atom physics, where atom-atom 
cross sections are calculated to three decimal places or 
better. The difficulty of calculating the exciton- exci- 
ton scattering cross section has nothing to do with de- 
tails of the band structures of the crystals, which are 
well understood. The primary difficulty is the fact that 
the electron and hole have comparable mass, so that all 
the following exchange effects are important and cannot 
be neglected: electron-electron exchange, hole-hole ex- 
change, and electron-hole exchange. In atoms, there is 
no electron-nucleus exchange, and nucleon-nucleon ex- 
change can typically be ignored (although M. Combescot 
has argued that under the conditions of BEC the nucleon- 
nucleon exchange may be important, since the wave- 
length, by definition, is comparable to the distance be- 
tween the particles; see Ref. 44] and references therein 
for the general approach of accounting for these exchange 
effects). 

The most accurate calculations for the exciton-exciton 
scattering cross section in CU2O have been done by 
Shumway and Ceperley [isj . and give a cross section of 
5 A, which is comparable to the Bohr radius of the ex- 
citons. In that work, the spin of the electrons and holes 
was assumed to be aligned. The spin structure of the ex- 
citons can lead to significant differences in the scattering 
cross sections between different spin states 46] , although 
these calculations have not been worked out with the 
same accuracy as the Ceperley and Shumway work. 



II. EARLY EXPERIMENTS ON EXCITON BEC 
IN CU 2 

A. Spectral analysis experiments 

The first experiments to show promise for exciton BEC 
in CU2O used excitation with a green laser which was 
absorbed within 5 /jm of the surface [UEs!]. This created 
a thin "pancake" of excitons which were not trapped, and 
could flow away from the surface, but on time scales of 
a few nanoseconds remained at high density at the point 
of creation. 

The phonon-assisted luminescence spectrum of the or- 
thoexcitons (the upper, triplet state of the Is excitons in 
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FIG. 2: Phonon-assisted luminescence spectrum showing 
the spatially-averaged energy distribution of orthoexcitons in 
CU2O for two different times as the density falls after the ex- 
citons are created by a short laser pulse. The fit temperature 
at 5 ns is approximately 60 K while the temperature at 15 ns 
is approximately 30 K. From Ref. [48| . 



C112O) could be well fit by the Bose-Einstein distribution 
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with a = l/i/ksTl ~ 0.1 over a wide range of densities. 
In other words, the gas remained near the phase bound- 
ary for BEC but did not cross it. This was termed "Bose 
saturation." Figure [2] shows that this saturation corre- 
sponds to an energy scaling such that the shape of the 
photoluminescence remains the same. This scaling was 
found to apply over more than an order of magnitude of 
density variation. 

The fits to the Bose-Einstein distribution imply an ab- 
solute density of the particles according to 



N = 



1 



o(E-/j.)/k B T _ l 



D(E)dE. 



(2) 



The densities from the fits were compared to data for the 
relative change of the density, found by taking the total 
luminescence intensity and dividing the measured exci- 
ton cloud volume at each point in time, and the variation 
in the density deduced from the fits agreed with the vari- 
ation in density from the intensity data within a factor 
of two over a wide range of density [I?} . 

Later experiments looked at the paraexciton phonon- 
assisted luminescence under very similar experimental 
conditions and concluded that the paraexciton density 



exceeded the density for BEC predicted by the ideal 
gas theory 49] . Line shape analysis of the paraexciton 
phonon-assisted emission was harder, because the emis- 
sion intensity of this line is 500 times weaker than the 
orthoexciton line, and it lies near another, brighter emis- 
sion line. The relative intensity of the paraexciton line 
seemed consistent with the conclusion that the paraexci- 
tons exceeded the critical density. 

The first indication that this analysis could not be 
the whole story came from noting that the many-body 
collision broadening of the lines (discussed above) was 
very low, and not consistent with the densities implied 
by the BEC fits [43j . Although the exciton-exciton scat- 
tering cross section is not known to better than a fac- 
tor of two, as discussed above, if the experiments were 
truly to have densities in the range of 10 19 cm -3 , as im- 
plied by the fits of Eq. (JTJ to data like that shown in 
Fig. along with the observed spectral broadening of 
less than 0.1 meV, then these two observations together 
would require a scattering cross section with radius much 
less than the Bohr radius of the excitons. The crucial 
evidence came in comparing the exciton luminescence 
spectrum in the case of a cold (2 K) lattice, when the 
exciton distribution fit a Bose-Einstein distribution at 
70 K, to the luminescence spectrum with_a 70 K lat- 
tice under the same excitation conditions 50] . The total 
number of photons emitted in each case was nearly the 
same, indicating that the densities were comparable, but 
the shapes of the spectra were quite different. The de- 
viation from a Maxwell-Boltzmann distribution in the 
cold-lattice case could therefore not be from Bose statis- 
tics at high density if the same density gave a Maxwell- 
Boltzmann distribution when the system was closer to 
equilibrium with the lattice. This pointed to spatial in- 
homogeneity in a nonequilibrium system as the cause of 
the non-Maxwellian distribution. 

The main ambiguity of these experiments comes from 
the fact that the luminescence from the exciton gas must 
be integrated over at least one dimension, since the pho- 
tons are collected by imaging a three- dimensional crystal 
in two dimensions. To fit the Bose-Einstein distribution 
to the data, some model of the spatial variation of the 
density and temperature must be used. This was recog- 
nized in the early work [481 ] , and it was found that nearly 
all models which assumed that the temperature and den- 
sity of the orthoexcitons decreased monotonically in the 
direction away from the surface gave spatially integrated 
spectra that looked very similar to the simple homoge- 
neous spectrum of Eq. (JTJ . (The difficulty of interpreting 
the three- dimensional data integrated over one dimen- 
sion was one motivation for pursuing experiments with 
excitons, or exciton-polaritons, in two dimensions (for 
Refs. [U and H3-) 

found that a reasonable fit to 
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Later modeling 



the spatially-integrated orthoexciton luminescence spec- 
tra could be obtained by assuming a large cloud of cold 
orthoexcitons far from the surface (see Fig. A sim- 
ple analysis would seem to indicate that orthoexcitons 
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do not travel very far within their lifetime due to down- 
conversion into paraexcitons [5^ - l56j . and therefore one 
should not expect so many orthoexcitons away from the 
surface. However, the model of Ref. [53| showed that 
picking a high enough rate for an exciton Auger recom- 
bination process could generate orthoexcitons at long 
distances from the surface due to up-conversion from 
paraexcitons. 

This brings up the Auger recombination process in 
CU2O, which was established to occur in CU2O in the 
1980's [36| and became an increasingly important topic 
of study in the 1990's. The Auger process occurs when 
two excitons collide, and one of them recombines, but in- 
stead of emitting a photon, the energy of the recombining 
exciton is given to ionizing the second exciton. The hot 
electron and hole thus produced can then lose energy by 
phonon emission or by collisions with other carriers, and 
finally form into an exciton again. It is assumed that 
the spins of the electron and hole are randomized in this 
process, so that the returning exciton can be either an or- 
thoexciton or paraexciton. The Auger process therefore 
1) reduces the lifetime of the excitons dramatically at 
high density, from the radiative lifetime of milliseconds 
down to less than a nanosecond, 2) heats the exciton 
gas due to the gap energy given to the ionized electron 
and hole, which then interact with the rest of the sys- 
tem, and 3) gives a population of orthoexcitons in excess 
of that expected from the Boltzmann occupation factor 
e -A/k B T ^ wnere a is the paraexciton- orthoexciton split- 
ting. Such an imbalance implies that the system is not in 
full equilibrium, which can occur when the lifetime due 
to Auger annihilation of the paraexcitons is shorter than 
the ortho-to-para conversion time. 

Throughout the 1990's and early 2000's, several groups 
tried to measure the rate of Auger recombination in 
CU2O. This work is reviewed below. The measurements 
are made difficult by the fact that there is also a two-body 
conversion process by which two orthoexcitons can collide 
and turn into two paraexcitons. Since this process has 
the same density dependence as the Auger process it can 
be difficult to distinguish the two processes in many ex- 
periments. In recent experiments, the group of Kuwata- 
Gonokami used infrared absorption measurements to ex- 
cite the ls-2p exciton transition [5Ti - f59j | . which is analo- 
gous to the ls-2p transition in alkali atoms. Since this 
transition only occurs when the excitons exist, the ab- 
solute absorption cross section of the transition gives an 
absolute measurement of the number of both species of 
excitons. Therefore the spin flip process, which conserves 
total number of excitons, could be distinguished from the 
Auger process, which annihilates excitons. 

Although there was a range of values found for the 
Auger rate, the general consensus arising from the Auger 
measurements was that the density of the excitons could 
not be as high as predicted by the fits of the surface- 
excitation data to the Bose-Einstein distribution, be- 
cause the lifetime would be too short at those densities. 
Since, as discussed above, there was also inconsistency 
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FIG. 3: Fits of a theoretical model for nonequilibrium flow of 
excitons to the phonon-assisted photoluminescence spectrum 
of CU2O with high-intensity surface excitation for three exci- 
tation powers under conditions similar to those of Fig. 
two different times after a short laser pulse. From Ref. 



in the broadening of the spectral lines and the abso- 
lute number of photons coming from the phonon-assisted 
luminescence, it was generally concluded by the end of 
the 1990s that the early experiments did not show Bose- 
Einstein condensation. The Bose saturation effect, which 
corresponds to an invariance of the luminescence spec- 
trum as temperature is scaled (see Fig. is therefore 
surprising, because there are a number of details which 
go into the classical numerical model of which have to be 
fortuitously related to give the nearly invariant spectrum 
which is observed; without these coincidences the spec- 
trum would have different shapes at different times. The 
model of Ref. [53| used choices of the physical parameters 
which were reasonable based on the experimental data, 
but even with the optimal choice of realistic parameters, 
the fit of the model to the data is qualitatively similar in 
shape but not especially good (see Fig. \S§ . The satura- 
tion effect corresponds to a scaling law, which we discuss 
below, in Section UlIBI 



B. Transport experiments 

Also during the 1990's, another set of experiments 
was performed that focused on transport of the excitons 
rather than spectral signatures [60j-|63j . The experiments 
were done as follows: a long-wavelength (red) laser with 
photon energy nearly resonant with the exciton energy 
was used to flood a whole three- dimensional crystal of 
CU2O with a low-density exciton gas. Then a very in- 
tense pulse of green light was used to excite one surface 
of the crystal, creating a very dense exciton cloud at the 
surface. Excitons were then detected on the opposite side 
of the crystal, in some cases over a millimeter distant, by 
a bimetal detector that ionized the excitons into free elec- 
trons and holes which then produced an external current. 
The signal would only occur if excitons migrated across 
the crystal into the detector. 
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The experiments showed that the exciton signal on the 
back side of the crystal jumped up in a sharp pulse when 
the green laser pulse intensity exceeded a critical thresh- 
old. The arrival time of the pulse corresponded to propa- 
gation at the speed of sound across the crystal. This was 
interpreted by the authors as superfluid motion of the 
excitons across the crystal when the excitons exceeded 
a critical density threshold. The high density pulse of 
excitons created by the green pulse was viewed as being 
amplified by the large population of cold excitons created 
by the resonant laser. 

There were soon alternative interpretations of these 
experiments. Reference (64j pointed out that the width 
of the pulses was orders of magnitude longer than the 
expected coherence length, and that a classical soliton 
model could reproduce many of the effects. Tikhodeev 
and coworkers |6a - [6a | pointed out that the intense green 
laser pulse also created many hot phonons which could 
blow as a nonequiilibrium "phonon wind" away from the 
surface, pushing excitons at the speed of sound, and rea- 
sonable models of the phonon wind could explain the ex- 
citon propagation. The phonon wind was well established 
(e.g., Ref. mM) 

in the 1980's as an effect in surface- 
excited solids which can push excitons and electron-hole 
liquid for hundreds of microns. One experimental result 
in support of this interpretation is that a second pulse was 
seen at the detector at much later times, corresponding 
to three passes across the crystal at the speed of sound, 
i.e., hitting the back surface, bouncing off, and coming 
back one round trip later [7l[ ■ This strongly supports the 
view that a sound pulse was created by the green laser 
pulse, and this sound pulse propagated across the crystal 
and reflected from the surfaces, sweeping excitons from 
the cold exciton gas as it traveled. 

As with the spectral signature experiments, it was 
widely concluded by the end of the 1990's that these 
experiments with surface excitation also did not show 
Bose-Einstein condensation. Both of these experiments 
used very intense surface excitation in order to get very 
high exciton density, but both had the same drawback 
that the exciton gas thus created was very far from equi- 
librium and unconstrained in space, so that flow out of 
the excitation region had to be modeled, including the 
heat flow and phonon wind. Because of this, attention 
returned to trapping the excitons using stress. 



C. A Way Forward 

It became clear by the first ICSCE conference in 2003 
[72| that the Auger effect played a dominant role in pre- 
venting EBEC in CU2O by shortening the lifetime of the 
excitons. To avoid this, the density of the excitons can be 
dropped, with a corresponding drop in the temperature, 
since T c of the BEC transition is proportional to the den- 
sity to the 2/3 power in three dimensions. At first, one 
might think that dropping the density would not help, 
because the collision rate between the particles has the 



same dependence on density as the Auger rate, and there- 
fore the ratio of the Auger lifetime to the collision time 
will not change. Excitons in a semiconductor have an- 
other means of thermalizing besides collisions, however, 
which is exciton-phonon scattering. The time scale for 
exciton-phonon equilibration is hundreds of picoseconds 
to nanoseconds [32|,[73j at low temperature, so that when 
the gas has low density, the time scale for thermalization 
via phonons can be much shorter than the lifetime due 
to Auger recombination. 

One issue in some solid-state systems in going to lower 
temperature and density is that the excitons can become 
localized in minima of the disorder (e.g., [74]])> so that 
they do not move and do not thermalize. Many CU2O 
samples have extremely good quality, however, so that 
disorder is not a limiting factor. Another issue is that the 
polariton effect, i.e., mixing of the exciton and photon 
states, becomes important for excitons at low momen- 
tum, i.e., low temperature. There is a polariton effect 
in CU2O for both the orthoexciton species and for the 
paraexciton species in stressed crystals; although the di- 
rect recombination process for paraexcitons is forbidden 
by symmetry in unstressed crystals, if stress is used to 
create a trap for the excitons, the paraexciton coupling to 
the photons becomes comparable to that of the orthoex- 
citons 36]. It is possible, however, to go to temperatures 
so low that the excitons have momentum which is below 
the mixing region of the photon and exciton states. 

The above considerations imply that the best temper- 
atures to see BEC effects in CU2O are in the sub-Kelvin 
range of tens to hundreds of mK. This removes one of the 
original motivations for studying exciton condensates, 
which is to see high-temperature condensation; the fact 
that the excitons in CU2O are stable up to room temper- 
ature turns out not to help in getting them to condense 
at room temperature. One can still hope to see a true ex- 
citon BEC in a three-dimensional system, however. Two 
experimental groups have pursued EBEC experiments of 
excitons in traps in CU2O at milliKelvin temperatures in 
the past five years; these are reviewed in Section HVl 



III. REVIEW OF THE AUGER PROCESS 

As discussed above, one of the most important of the 
mechanisms which determine the lifetime of excitons in 
CU2O is believed to be a nonradiative two-body decay 
process, as evidenced by numerous experiments, which 
have shown that the exciton lifetime is density depen- 
dent. This "Auger" process is a well-known decay mech- 
anism in the physics of semiconductors. This process is 
also present in electron-hole plasmas, in which the recom- 
bination of an electron and a hole excites either an elec- 
tron high in the conduction band or a hole deep in the va- 
lence band. This process, which occurs in many semicon- 
ductors, has long been studied both theoretically [7oT - [78| 
and experimentally (79l |8C| . A similar decay process has 
also been observed in experiments with cold metastable 
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He atoms, where the term "Penning collisions" is used 
[8lj |. It also occurs with Frenkel excitons in organic ma- 
terials, where it can be quite strong [H, There may, 
in fact, be a general principle that excitons with large 
binding energy, and therefore are Frenkel-like, intrinsi- 
cally will have strong Auger recombination. 

In CU2O, the exact nature of the Auger process is still 
an open question and is under current investigation. The 
general picture is that in such a process two excitons 
collide, the one recombines, transferring its energy to the 
other, which ionizes. The electron and the hole which 
result from the ionization process then form an exciton, 
after a fraction of their energy is transferred to the lattice, 
via phonons. Therefore the Auger process is not only a 
decay process, but it heats the system as well, since the 
exciton energy is transferred partly to the exciton gas 
and partly to the lattice. 

Before we review the main results on this process, we 
start with some remarks and also mention some reasons 
which make the study of this process in CU2O difficult: 

(1) The experiments which have measured/examined 
this process cover a rather wide range of (exciton) tem- 
peratures T (which is not necessarily equal to the lattice 
temperature) and exciton densities n. Thus, in trying to 
extract the temperature and density dependence of this 
process, it could well be that its rate has a different func- 
tional dependence on T and n, depending on their actual 
value /range. 

(2) In excitons in CU2O there are many different pro- 
cesses which take place at the same time (radiative re- 
combination, conversion between the orthoexcitons and 
the paraexcitons, expansion, etc.) These processes ob- 
scure the measurement of the Auger decay process alone, 
since there is strong evidence that the rate of some of the 
above mechanisms is comparable to the rate of the Auger 
process. 

In particular, the interconversion process of the to- 
tal angular-momentum triplet state orthoexcitons into 
the total angular-momentum singlet state paraexcitons 
is of great importance. It is believed that it takes place 
via two mechanisms, namely a phonon-assisted process 
54l [56l [Hj-86] and a collisional spin-exchange process 
87j, as we discuss below. At low temperature and high 
density, the collisional ortho-para conversion process in 
many ways mimics the Auger collisional recombination 
process, making analysis much more difficult. 

As a result of the many processes which all take place 
at the same time, one has to write down rate equations, 
which then give the Auger rate (as well as the rate of the 
other processes) from some fitting. Therefore, to a large 
extent, the determination of the Auger rate is model- 
dependent, and this complicates the problem. In many 
cases, the models must include parameters for spatial 
diffusion of the excitons in order to estimate the volume 
of the exciton cloud. 

In some regimes of parameters, different processes can 
be isolated. For example, experiments can be done in a 
very low density regime in which collisions between ex- 



citons are negligible. In this case, the phonon- assisted 
conversion between orthoexcitons and paraexcitons can 
be determined very accurately [H], HU . It is also possi- 
ble to establish the intrinsic lifetime of the excitons by 
fitting the decay of the luminescence at very low densi- 
ties, when the Auger recombination rate is negligible. In 
addition, at temperatures which are high compared to 
the energy splitting between the ortho and para states, 
chemical equilibrium between the ortho and para popu- 
lations can be assumed, so that all ortho-para conversion 
processes drop out of the rate equations. This fact has 
been used to extract the total Auger rate of the exciton 
population in the high-temperature regime [H| , although 
these results still relied on a model for the exciton cloud 
volume. 

(3) The determination of the exciton density has 
proven to be a rather difficult task. There are basically 
three approaches that one may follow to extract the den- 
sity. The first is the spectroscopic approach, in which 
some specific phonon-assisted recombination line is fitted 
to a Bose-Einstein distribution and from that the density 
and the temperature are extracted [47J , as discussed in 
Section [1] above. As discussed above, the estimates of 
density based on spectroscopic fits to an equilibrium ho- 
mogeneous gas strongly overestimate the density. The 
second method of extracting the exciton density is based 
on estimating the volume of the exciton gas and the num- 
ber of excitons [13, [H, HH . The third method relies on 
the absorption spectrum that corresponds to the Is to 
2p radiative excitonic transition. This method does not 
depend on the strength of the recombination lines and is 
thus advantageous in that respect. In addition, as we dis- 
cuss in more detail below, the absorption spectrum of this 
process is very sensitive to the degeneracy of the excitons 
and it can serve as a clear and indisputable probe for the 
transition to a Bose-Einstein condensed phase [|cj ■ 

(4) The theoretically-calculated decay rate involves 
some parameters with relatively high uncertainty (in 
their numerical value) raised to some large powers. As a 
result, the theoretical estimate of the decay rate is also 
rather uncertain. 

The importance of the Auger process in excitons in 
CU2O has been revealed in several experiments over the 
last three decades [H HI H M, H, HI, HI M, M~ 
[IH . The inverse of the total decay time of excitons shows 
an approximate proportionality to their density, which is 
an indication of a two-body collision process. The Auger 
process also leads to the presence of orthoexcitons in 
stress-generated traps for times much longer than the or- 
thoexciton lifetime; although the exciton gas in the wells 
consists primarily of paraexcitons, high energy electrons 
and holes produced by the Auger ionization of para ex- 
citons lead to re-formation of excitons in essentially ran- 
dom internal angular momentum states, and thus to both 
ortho and paraexcitons [H, [93| . 

The difficulties discussed above have led to a wide scat- 
ter in experimental estimates of the Auger decay rate. 
The decay rate has been extracted to be negligible |)5[ 
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and also non- negligible, and when non-negligible, to be 
independent of the temperatur e 1591 . to increase linearly 
with the exciton temperature (88], and to be inversely 
proportional to the exciton temperature (86j . The stress 
dependence of the Auger decay rate is also of impor- 
tance [55j. This process is sensitive to the various sym- 
metries of the crystal, which in turn impose certain selec- 
tion rules/constraints on the matrix elements and on the 
decay rate, while stress tends to lift these constraints. It 
has also been argued [%[ that excitons form biexcitons, 
which then rapidly decay via the Auger process, however 
there is no direct experimental evidence for this process, 
at least to our knowledge. In addition, the theoretical 
calculation presented in the same reference ignores the 
orthogonality of the bands of the crystal (according to the 
authors of this study the Auger mechanism is thus associ- 
ated with breaking the band symmetries, e.g. by impuri- 
ties, and therefore differences of impurity concentrations 
could explain the different Auger rates measured). 



A. Review of Auger theory, constrained by 
experiments 

Turning to the theoretical study of the Auger decay 
rate, according to Ref. (9?J this process may be either "di- 
rect" , or phonon assisted. In the phonon-assisted process 
there is a phonon involved, which is absent in the direct 
one. While the participation of a phonon suppresses the 
rate of this process, on the other hand, it lifts the sym- 
metry constraints set by the band structure of CU2O. 
We should stress here that there is an intimate connec- 
tion between the radiative and the Auger decay mech- 
anisms: the radiative lifetime of excitons is rather long 
as the dipole matrix element between the conduction and 
the valence bands (which have the same parity) vanishes. 
Phonon-assisted mechanisms provide an alternative way 
for the excitons to decay radiatively. As in the Auger pro- 
cess, the exciton-phonon interaction suppresses the rate 
of this process, but on the other hand it makes the process 
dipole allowed. As a result, the phonon-assisted processes 
are actually the dominant radiative decay mechanisms. 
Rather similar results also hold for the direct and the 
phonon-assisted Auger processes, as we discuss below. 

The excitons in CU2O in the condensation experiments 
are formed from electrons in the lowest conduction band 
(rjj") and holes in the highest valence band (F^) (the 
"yellow series" ) . The dependence of the matrix element 
for the Auger collisions on the momentum exchange be- 
tween the colliding excitons is of major importance [97| . 
If this is constant, then the rate is proportional to the 
density for low temperatures and its only temperature 
dependence comes from phonon emission [88| . If the ma- 
trix element is proportional to the momentum exchange, 
the Auger decay rate is proportional to the exciton tem- 
perature and to the density of the non-condensed exci- 
tons. 



fluenced by the following factors: (i) the same (positive) 
parity of the lowest conduction and the highest valence 
bands of CU2O, which form the yellow excitons (ii) the 
fact that the valence band is not a pure spin state and 
(iii) the symmetry and the location of the negative- parity 
bands with respect to the conduction and valence bands 
which form t he y ellow excitons. According to theoretical 
calculations [88|, |97[ the dominant Auger process is ac- 
tually the phonon-mediated process, very much like the 
radiative recombination process. As we mentioned also 
above, although in this case one also needs a phonon, on 
the other hand, there is no symmetry constraint, which 
enhances the rate of this process. 

Based on Ref. [I?]], Ref. [88[ has evaluated the Auger 
rate T^ uger — An D for orthoexciton-orthoexciton colli- 
sions, where 



A w 3.698 x 1CT 17 ( 1 
+1.298 x 1CT 17 ( 1 



e 160/T _ l 
2 

e 217 / T _ 1 



(3) 



According to Ref. |97J , the Auger process is strongly in- 



Here the rate T^ r is measured in nsec -1 , the tempera- 
ture T is measured in Kelvin, and the orthoexciton den- 
sity n in units of 10 16 cm -3 . The two exponential factors 
come from the emission of the optical phonons T^ 2 and 
IYg, with energies 13.8 meV and 18.7 meV, respectively. 
The above formula is in good agreement with the experi- 
mental measurement of this process [88[ (where the den- 
sity was determined spectroscopically) for a wide range 
of temperature between 2 K and 235 K with essentially 
one fitting parameter, which is reasonably close to its 
expected value. According to Eq. ([3]), for T roughly up 
to 77 K the rate is independent of T, while for higher 
temperatures the rate increases linearly with T (which is 
obvious if one expands the exponentials) . 

References [89| and [HJ have estimated the exciton 
number and the volume of the gas at relatively low tem- 
peratures (and thus did not get the density spectroscop- 
ically). The value of the Auger constant A that they 
found using this different method was 7 x 10 -17 cm 3 /ns 
and 6 x 10~ 17 cm 3 /ns, respectively, which is in decent 
agreement with the result of Eq. ([3]) and with the exper- 
imental results of Ref. (88|. Reference 31 1 also recently 
measured the orthoexciton-orthoexciton Auger rate to be 
5 x 10~ 17 cm 3 /ns at low temperature, very close to the 
low-temperature value of Ref. [H| . 

For temperatures above 100 K, it is worth a close 
look to compare the results of Refs. [86[ and [88J. As 
noted above, the values for the Auger constant of the two 
groups at low temperature fall within the same range, on 
the order of 10~ 16 cm 3 /ns, with about a 30% deviation, 
which is not surprising given the various experimental 
uncertainties. However, in the range 100-200 K, the two 
experiments showed opposite trends with temperature: 
the Auger rate was found to rise by about a factor of 
two in Ref. [88[ and to fall by a little less than a factor 
of two in Ref. 86]. The two experiments used similar 
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methods, with pulsed near-resonant dye laser excitation. 
Even apart from any modeling, the raw data of Ref. [88j 
clearly showed that the excitons had shorter lifetime at 
high temperature, while the raw data of Ref. [§6[ shows 
the opposite trend. 

The primary difference between the experiments was 
that a much tighter focus was used in the experiments of 
Ref. |86|, about 15 /zm radius, as compared to 120 /im 
radius in Ref. |88j |. This gave higher densities and much 
greater range of lifetime, but may also have had other ef- 
fects. One possibility is that the local lattice temperature 
may have gone well above the background bath tempera- 
ture due to the concentration of energy, leading to a local 
"hot spot" (cf., e.g., Refs. [Hill.) The exciton temper- 
ature may also have gone well above the lattice tempera- 
ture, if the Auger rate was fast compared to the phonon 
emission rate. Taking the Auger rate as increasing with 
temperature as discussed above, if a local exciton cloud 
is hotter than the rest of the crystal, a higher Auger rate 
will occur. To actually have a decrease of the Auger rate 
with increasing lattice temperature, however, would re- 
quire that the local temperature excursion was greater at 
low bath temperatures. It is also possible that in the ex- 
periments of Ref. [88| , at the highest temperatures there 
was also an excursion of the temperature above the lat- 
tice temperature, because the absorption length of the 
input laser shortened to approach 20 /xm (although the 
lateral size still gave much lar ger volume for the spot size 
than the experiments of Ref. [86|). 

Alternatively, or possibly in conjunction with a local 
hot s pot , the rapid expansion of the exciton cloud seen in 
Ref. 86| may have played a role. In the experiments of 
Ref. [88j . volume expansion was a minor factor because 
the initial size of the cloud was large. With the smaller 
initial spot size in Ref. [86| . the expansion of the cloud 
played a bigger role. The authors of Ref. [8(| carefully 
estimated the volume as a function of time V(t) based on 
a model of diffusion to get the density n(t) = N(t)/V(t), 
but the relative uncertainty of the volume will be larger 
with a smaller spot size. 

Atpresent, Refs. [5£| and [9j| are the outliers, with 
Ref. [59j reporting a much larger value for the Auger rate 
than the rest of the literature, approximately 4 x 10~ 16 at 
low temperature, and Ref. [95] reporting a much lower, 
indeed negligible, rate, of the order of 10 -23 cm 3 /ns. The 
majority of other experimental work [3ll l86l . [88j j gives 
Auger rates in the range of 10~ 17 to 10~ 16 cm 3 /ns. 

In the case of the work by the group of Gonokami, 
a multi-step process of analysis and measurements was 
used to derive the Auger rate. In the first step, the ab- 
solute number of excitons is established by far infrared 
Lyman absorption spectroscopy of the ls-2p, ls-3p, etc. 
transitions for both orthoexcitons and paraexcitons [571 - 
[59[ . The reason that the Lyman spectroscopy can estab- 
lish the absolute number of excitons while band-to-band 
absorption cannot is that the absolute cross section for 
the Lyman series can be very accurately calculated from 
first principles for intraband transitions, while the band- 



to-band oscillator strength cannot. 

Having this method to establish the absolute number 
of both ortho and paraexcitons, the number must be con- 
verted to a density by an estimate of the volume of the 
exciton cloud. In a stress-generated harmonic trap, this 
can be done by imaging the exciton cloud from two sides 
of the crystal, but typically the number of excitons in 
such a trap is too few to generate a reasonable Lyman 
absorption measurement. Therefore a large exciton cloud 
was generated in a bulk crystal, without a trap, and a 
model for the exciton diffusion was used to estimate their 
volume. This modeling is the greatest source of uncer- 
tainty in the measurements, since various effects such as 
phonon wind, discussed above, may also affect the exci- 
ton motion. The authors of Ref. [59| used the d iffusion 
constant for paraexcitons reported in Ref. |l00j |. which 
was very high. This was reasonable, given that their 
sample was cut f rom the same natural sample as that 
used in Ref. [l00j |. but was not checked directly. If the 
diffusion constant was estimated too high, then the vol- 
ume would be estimated as much greater, and therefore 
the density much lower, which could give systematically 
high Auger rates. 

The work of Ref. 59] found an Auger rate which 
was inde pendent of temperature below around 70 K. In 
Ref. [l09j | , the authors argued based on this that the cross 
section for Auger recombination must diverge at low tem- 
perature proportional to 1/v, where v is the velocity of 
the particles, since in the semiclassical Born approxima- 
tion in three dimensions, a collisional process has rate 
1/t = nav, where n is the density and a is the cross 
section. It is not clear, however, that the Born approx- 
imation and the consequent concept of cross section ap- 
plies here. As discussed above, Eq. ([3]) predicts an Auger 
recombination rate independent of temperature below 
around 70 K, which should be the same whether the ex- 
citon gas is condensed or not. 

The work of Ref. j9f| gave the opposite result of very 
low Auger rate for both paraexcitons and orthoexcitons. 
In this experiment, the absolute density of the excitons 
was estimated by using the effect of screenin g of e xcitons 
at high density on the absorption spectrum 10l|, based 
on a theory of the Mott transition. Their theory of the 
screening led them to estimate initial densities on the or- 
der of 10 19 cm -3 , and under these conditions they saw 
no evidence of an Auger effect. The theory of the exci- 
tonic Mott t ransition is hotly debated and quite complex, 
however [l3|, Il02l Il03| , and therefore it is difficult to find 
agreement on a universal theory for the density at which 
the exciton absorption will be screened out. 

The above experiments primarily were sensitive to the 
ortho-ortho Auger recombination rate. The question of 
the value of the Auger collisions between paraexcitons 
is even more more subtle. Reference [3l| has provided 
evidence for a very low paraexciton-paraexciton Auger 
rate, on the order of 10 -18 cm 3 /ns. Earlier experiments 
of Refs. (5f| in stressed crystals also gave a para-para 
Auger rate about two orders of magnitude less than the 
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ortho-ortho rate. However, other experiments suggest 
that Aug er collisions are also possible among paraexci- 
tons [89f . Reference [HI] used a single Auger rate con- 
stant for all collisions in their model, so that their fit 
was not sensitive to a difference between the para-para, 
ortho-para, and ortho-ortho rates. 

Theoretically the para-para Auger collisions should be 
negligible [97[ at zero stress (the effect of stress is dis- 
cussed below). The reason for this is the band structure 
of CU2O, which makes the direct Auger process for para- 
para collisions forbidden, while the phonon-assisted pro- 
cess is highly suppressed [97] . This is due to the fact that 
the phonon-assisted Auger process for para-para Auger 
collisions goes through a deep valence band, as opposed 
to the ortho-ortho collisions, for which this process goes 
though a conduction band which is close in energy to the 
lowest conduction band. 



B. Orthoexciton to paraexciton conversion process 
and "saturation" of the orthoexciton gas 
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As discussed above, at very low densities, when the 
two-body processes are negligible, it has been observed 
that orthoexcitons convert into paraexcitons through 
a phonon-assisted process [84 - [86j . According to 
Ref. [54| its rate scales as T 3 / 2 . More recently Refs. 85 ] 
and [86[ have reported that its rate is constant at very 
low temperatures and increases linearly with the temper- 
ature at higher temperatures. 

In addition to the phonon-assisted process, orthoexci- 
tons convert into paraexcitons via a spin-exchange pro- 
cess, where two orthoexcitons collide, they exchange their 
electrons or holes, resulting into two paraexcitons in the 
final state. This process was studied theoretically in 
Ref. [13] , while Ref. [HI has shown experimental evidence 
for it. The decay rate of this process is proportional to 
the orthoexciton density n Q . 

A simple argument has been proposed at v arious times 
over the years (see, e.g., Refs. [l04j and jl05| ) to explain 
the saturation effect of orthoexcitons, discussed in Sec- 
tion [U which corresponds to a power law n oc T 3 / 2 . 
Recall that the exciton gas, which in general has a 
higher temperature than the lattice temperature, ex- 
changes phonons with the lattice and thus loses energy, 
while on the same time it is heated by the Auger process 
and the ortho-to-para conversion. From simple deforma- 
tion potential theory, the rate of loss of energy is pro- 
portional to T 3 / 2 (see Ref. 12], Section 5.1.4). On the 
other hand, both the Auger process and the collisional 
conversion process of orthoexcitons to paraexcitons heat 
the gas. In both cases the rate of energy gain of the gas is 
bn , where n D is the orthoexciton density. As a result, we 
can write a rate balance equation for the entropy of the 
orthoexcitons divided by the orthoexciton number S /N 



dt \N 



-aT 3/2 + bn , 



(4) 



FIG. 4: Energy loss rate for orthoexcitons in CU2O as a func- 
tion of exciton temperature, for a cold lattice. Dashed line: 
rate when only acoustic phonon emission is included. Dotted 
dashed line: rate when acoustic-phonon and optical single- 
phonon emission are included. Solid line: rate when acoustic, 
optical single-phonon, and optical two-phonon emission are 
included. Dotted line: the T 3 ^ 2 power law. From Ref. [3^ ]. 



which implies that the orthoexcitons move along adia- 
bats, where n oc T 3 / 2 , in agreement with the experi- 
mental data [93l | . This above argument a pplie s both in 
the classical and in the degenerate regime [l05j ] . 

A problem with this argument, however, is that the 
rate of heat loss for orthoexcitons in CU2O never exactly 
follows the T 3 / 2 law predicted by elementary deforma- 
tion potential theory. There are both high-temperature 
corrections, due to optical phonon emission and high- fre- 
quency corrections of the exciton wave function, and low- 
temperatur e cor rections due to momentum-conservation 
limitations [lOOj ]. that give deviations from this law. The 
rate of phonon emission in CU2O has been measured care- 
fully and reported in Ref. [32j ; a summary plot is shown 
in Fig. |U 



C. The effect of stress 

As we mentioned earlier, the Auger decay process is 
strongly influenced by stress. In stressed crystals the 
selection rules which determine the Auger decay rate are 
no longer valid, and as a result stress enhances the rate 
of this process. 

The most important issue in the case of stressed crys- 
tals is the effect of stress on the paraexciton-paraexciton 
Auger collisions. As discussed above, in non-stressed 
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crystals this is an open question, however it is certain that 
under stress these collisions are allowed. An analogous 
situation is the radiative decay of paraexcitons, where 
stress makes the direct radiative recombination process 
of paraexcitons allowed, and it also enhances the phonon- 
assisted recombination lines. More specifically, uniaxial 
stress mixes the and T% valence bands with the re- 
sult that the para-exciton energy increases for small val- 
ues of the applied stress but decreases for higher values, 
and the direct radiative recombi nation inc reases quadrat- 
ically with increasing stress [HI, Il06l4l08| | . 

Indeed, the experiment in Ref. [36| shows clearly that 
para-para Auger collisions are allowed in stressed crys- 
tals. Uniaxial stress makes the direct Auger process al- 
lowed, and also the phonon-assisted processes. In both 
cases the orientation of the uniaxial stress is very im- 
portant, since it d etermines the selection rules in the de- 
formed crystal pi ll06Hl08| . These processes are respon- 
sible for the presence of orthoexcitons in the quantum 
well of Ref. [36|, that is created by the applied stress. 
The orthoexcitons that are present are generated by the 
reformation of Auger-ionized paraexcitons in random an- 
gular momentum states. 

The more recent experiment of Ref. [55[ has found 
that the Auger recombination rate for paraexciton col- 
lisions increases, while the rate of conversion of orthoex- 
citons into paraexcitons decreases. Furthermore, the 
same study has found that for zero stress the paraex- 
citon Auger collisions are negligible. As discussed above, 
these results are consistent with the theoretical study of 
this problem. 

The fact that the paraexciton Auger process becomes 
stronger under crystal stress creates a drawback for using 
stress to create a harmonic trap for the excitons. While 
the trap can be used to increase the density of the ex- 
citons, it also shortens the lifetime. Which process wins 
out will depend on the temperature, the trap depth, and 
other experimental details. 

D. Effect of the indistinguishability 

We mentioned earlier that according to theoretical cal- 
culations the paraexciton-paraexciton Auger collisions 
are expected to be negligible. Still, paraexcitons may 
experience Auger collisions due to the presence of or- 
thoexcitons. 

In addition, the expression of Eq. (J3J refers to the 
Auger collisions without including the effect of indistin- 
guishability in the orthoexciton-orthoexciton and in the 
paraexciton-paraexciton Auger collisions. According to 
Ref. (97[ > since the Coulomb interaction does not produce 
any spin- flips, 




where N a and N p are the numbers of ortho and para exci- 
tons, respectively. Recalling that in the Auger collisions 
one exciton is lost, while the ionized electron and hole 
form an exciton in random angular momentum states, 
the net loss of ortho and para excitons due to the Auger 
process is 
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IV. RECENT EXPERIMENTS ON EXCITON 
BEC IN CU a O 

Two currently ongoing experiments are attempting to 
produce EBEC in CU2O under the most optimal condi- 
tions. In these experiments, stress is used to produce a 
three-dimensional harmonic trap for paraexcitons. The 
lattice temperature is taken to the milliKelvin range, so 
that paraexcitons can equilibrate in states below the po- 
lariton mixing region, and also so that very low densities 
can be used, giving a very slow Auger recombination rate. 

As discussed in Section MI A[ the group of Gonokami 
in Japan performed measurements of the Auger recombi- 
nation rate of the excitons in CU2O using exciton Lyman 
spectroscopy to establish the number of excitons and a 
theoretical model for the diffusion to estimate the volume 
of the exciton cloud. Based on this method, the Auger 
cross section was found to be independent of temperature 
over the range of 4-70 K 59] . The authors suggested that 
this is similar to the case of inelastic collisions of atoms, 
but as discussed in Section III A, this result is also con- 
sistent with the phonon-assisted mechanism for Auger 
recombination presented in Ref. [971 ] . 

The authors then dropped the temperature into the 
range of tens to hundreds of milliKelvin and performed 
both spatial imaging and spectroscopy of the exciton pho- 
toluminescence. Initially, they found that the exciton 
cloud size agreed well with the predictions of equilib- 
rium of the exciton cloud with the lattice temperature. 
The cloud size depends on temperature in a trap because 
the harmonic potential of the trap leads to the relation 
r 2 oc fcsT. 

Be low a critical temperature, an "explosion" was found 
[l09j ] . in which the cloud expanded and was very far out 
of equilibrium (see Fig. [S]). One mechanism which could 
in p rinci ple lead to this is that proposed by Hijmans et 
al. [l lOl j for atomic hydrogen condensates. In this sce- 
nario, in a harmonic trap, the condensate may have such 
small volume and high density that density-dependent 
loss processes are strongly enhanced. The size of the 
condensate cloud is controlled by the repulsive interac- 
tions of the particles; if these are weak, then the size of 
the condensate can be quite small. 

Strange nonequilibrium effects of the paraexcitons at 
high density have been seen before, namely in the work 
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FIG. 5: "Explosion" of paraexcitons below a critical tempera- 
ture. The white dashed line indicates the trapping potential. 
The temperatures is fixed at 345 mK and the number of exci- 
tons in the trap was varied, estim ated at a) 2 x 10 7 , b) 5 x 10 s , 
and c) 2 x 10 9 . From Ref. [lol. 



done with surface excitation in the late 1980's and early 
1990's 48]. It is possible that when the paraexcitons 
approach the BEC state, some nonequilibrium recombi- 
nation effect does indeed prevent them from condensing. 

The German collaboratio n of Stolz and coworkers [3l[ 
and Frohlich and coworkers |l!3f also has seen equilibra- 
tion of paraexcitons in a trap, including at temperatures 
in the sub-Kelvin range. No evidence of an "explosion" 



of the type seen by Gonokami and coworkers was seen 
at low temperature, but the conditions were not exactly 
the same. In general, the degree of equilibration will 
depend on the rate of heat removal from the sample in 
the cooling system, the excess energy of the pump laser 
which creates the excitons, the intrinsic lifetime of the 
excitons, which depends on the impurities in the sample, 
and the duration of the laser pulse. Stolz's group has 
argued [lllL Ill2| that the mean-field renormalization of 
the exciton self-energy at high density should lead to a 
blue shift of the paraexciton line at high density and a 
change of the condensate spatial profile. Such a shift has 
not generally been seen in excitons in CU2O at high den- 
sity (see, e.g., Fig. H}. Although the mean-field density 
dependence is generally expected to lead to a blue shift 
(see, e.g., Ref. [12], chapter 8), and this has been seen 
with microcavity polariton condensates [114] and with 
indirect excitons in GaAs quantum double wells [ll5j j . 
in bulk semiconductors like CU2O the effect of higher- 
order correlations seems to cancel out most of this blue 
shift. A blue shift may indicate that much higher effec- 
tive densities have been reached, although also at higher 
temperature than needed for condensation. 



V. FUTURE PROSPECTS 

The process of the "explosion" seen by the Gonokami 
group is not well understood. It is possible that fur- 
ther theoretical study will allow a full understanding of 
this effect and its connection to condensation. In the 
meantime, the Gonokami group has reported that the 
explosion effect seems to be less dramatic at lower tem- 
peratures, and so experiments to go to even lower tem- 
peratures, which imply lower densities for condensation, 
may overcome this barrier. 

In a possible experiment that has been proposed and 
has been studied theoretically in Ref. [901 ]. it was shown 
that the absorption spectrum of infrared radiation induc- 
ing internal transitions of the excitons from the Is to the 
2p level is very sensitive to the degree of quantum degen- 
eracy of the gas. Actually, a similar idea has been used 
in the experiments with atomic hydrogen |l!6j ; however, 
in this case the mass of the excitons in the Is and the 2 p 
states is different (due to band structure effects |l!7| ). 
as opposed to the case of hydrogen, where the mass in 
the Is and in the 2p states is the same to a very good 
approximation. Remarkably, such experiments have been 
performed, as described above; see, e.g., Rcfs. [58, 59, 95]. 

In such an experiment one should observe the contri- 
bution of the orthoexcitons and of the paraexcitons to 
the absorption separately, with an energy separation of 
order A, assuming that the width of each distribution is 
of the order of k B T < A). 

As shown in Ref. [90j, the appearance of two distinct 
peaks in the absorption spectrum of ortho/para exci- 
tons would signal the presence of a Bose-Einstein con- 
densate, since in the condensed phase one deals with a 
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two-component system, and the two peaks would indi- 
cate the two different collective modes of it. Even if the 
ortho/para excitons have not crossed the phase bound- 
ary, but are highly degenerate, this would still show up 
clearly in the absorption spectrum. An advantage of this 
method is that it does not depend on the strength of the 
radiative recombination lines (which is very weak for the 
paraexcitons, as we mentioned above). 

As noted above, in practice, getting reasonable infrared 
absorption measurements of the ls-2p transition for a 
small number of excitons in a trap is difficult. It is pos- 
sible, however, to imagine using a stimulated process to 



enhance the absorption. 

Despite many false trails, research on excitons in CU2O 
has progressed and is still active, and there is good reason 
to believe that progress on establishing BEC of excitons 
in this well-studied crystal may still happen. 
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